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We show that low-energy photon scattering from a body with radius R slightly larger than its 
Schwarzschild radius resembles black-hole absorption. This absorption occurs via capture to one 
of the many long-lived, densely packed resonances that populate the continuum. The lifetimes and 
density of these meta-stable states tend to infinity in the limit r, — R. We determine the energy 
averaged cross-section for particle capture into these resonances and show that it is equal to the 
absorption cross-section for a Schwarzschild black hole. Thus, a non-singular static metric may trap 
photons for arbitrarily long times, making it appear completely 'black' before the actual formation 
of a black hole. 
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I. INTRODUCTION 

In this work, we consider the scattering of photons by 
the gravitational field of a non-rotating, finite-sized body. 
Such finite-sized bodies have radius R that slightly ex- 
ceeds their Schwarzschild radius Vg = 2GM/(?. We show 
that low-energy photon scattering from such objects re- 
sembles black hole absorption. The absorption arises due 
to the existence of a dense spectrum of narrow resonances 
(meta-stable states) whose lifetime and density tend to 
infinity in the limit Ts — > i?. 

Photons captured to such a resonant state are trapped 
on the interior of the body for a time t ^ h/Tn, where r„ 
is the width of a given resonance. For — >■ R, both the 
energy spacing D and width r„ tend to zero {t — )• oo), 
while their ratio remains finite. This allows us to de- 
fine the total cross-section for particle capture into these 
long-lived states using the optical model [l| , which is cal- 
culated by averaging over a small energy interval con- 
taining many resonances. At low energy the resonance 
capture cross-section is CTq = 47re^r^/3 where e is the 
energy of the incident particle. The absence of a longi- 
tudinal mode for photons means there is no state with 
total angular momentum equal to zero, hence the above 
cross-section tending to zero for zero energy. 

This cross-section exactly matches the absorption 
cross-section for a Schwarzschild black hole calculated by 
assuming complete absorption at the event horizon (see 
Refs. 12.-^1 ) derived previously for massless spin-1 parti- 
cles Note, however, that our calculation does not 
impose any special conditions at the boundary. Thus by 
considering the purely elastic cross-section of low-energy 
incident photons, we find that the absorption properties 
of a body with R resemble those of a black hole. 

The gravitational field of these near-black-hole objects 
is described using a suitable metric to model the interior, 
which is then joined to the standard Schwarzschild exte- 
rior metric at the boundary of the body. The above result 
is shown to be valid for any interior metric satisfying the 
following conditions: continuity with the Schwarzschild 
exterior at r = i?, a potential that deepens and develops 
a singularity in the black hole limit (allowing the particle 



to be treated semi-classically) and that this potential can 
be modeled as a harmonic oscillator in some finite region 
around the origin. We expect these conditions to be met 
by a large class of metrics; in this paper we present two 
metrics for which these conditions hold. 

Additionally, the results presented include the spin- 
(scalar) case considered previously (for j = 0) 
but given here for arbitrary angular momentum j. In 
this case taking Ts — >■ i? once again yields the black hole 
absorption cross-section given in Q . 

For a more detailed discussion of previous calculations 
involving scattering of various particles in the black hole 
spacetime we turn the readers attention to the introduc- 
tions in 0, [m and references therein. As in previous 
work fn\ we perform both analytical and numerical cal- 
culations, with good agreement between the two. 

II. WAVE EQUATIONS 

The Klein-Gordon equation for a massless spin-0 par- 
ticle on a curved manifold (with h — c = 1) reads 

d^iV^g^-d.^) - 0. (1) 

For massless spin-1 particles (photons) Maxwell's equa- 
tions on a curved manifold free of charges read 

dp{V^F'^0)^O, (2) 

where F"^ = d^A^ — and Af^ is the contravariant 

electromagnetic 4-potential. One may write a general 
static, spherically symmetric metric in the form 

ds^ = -e^'^-Ue + e^('-)dr2 + r^dn\ (3) 

By applying the separation of variables 

A, = e-^r^ir)sin9'-^p^ 

da 

to (m and ([2|) respectively in the metric ([3]), one derives 
at the radial differential equations for spin-0 and spin-1 
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respectively. In the latter we have implemented a gauge 
such that A^p is the only non-zero component [12i |. 

Defining s as the particle spin, either or 1, the radial 
equations for both spins can be combined to yield 
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v' -X' 
2r 



d(j) 
dr 



j(i + l)e 



(r) = 0, (4) 



where j the total angular momentum such that j — s >Q. 



III. GENERAL INTERIOR SOLUTION 

Equation (|4]) can be transformed into a Schrodinger- 
like equation by making the substitution (^(r) = x(r)/r 
and then mapping the radial coordinate to the Regge- 
Wheeler "tortoise" coordinate r* defined by dr* = 
e(A-'')/2^j,_ rpj^jg gives us 



1) d 



,u{r)-\(r) 



2r dr 



X(r*) = 0. (5) 



For a near-black-hole interior metric, e"''') for 
< r < i?, as time slows down in the limit — )■ R. 
This means that we may neglect all but the first term 
in square brackets in Eqn. ([5]) for all except very small 
energies. The solution regular at the origin is then 
x(r*) ~ sin(er*). For rs ^ R there is still a non-zero con- 
tribution to Eqn. dH) from the potentially large j{j + 1) 
term. This can be accounted for by introducing 93 as a 
constant phase shift determined by the boundary condi- 
tions at the origin (see below), giving 
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(6) 



which is now no longer regular at the origin, but is valid 
far from it. 

For small r, the j{j + 1) term in Eqn. (|4]) dominates 
and gives rise to the standard centrifugal barrier term 
assuming that, for small r, e'^ « const and e'^ — > 1. 
These assumptions are valid for the specific metrics we 
consider in Sec IIXI and apply to a larger class of static 
solutions where the potential is harmonic near the origin. 

The phase shift induced by the centrifugal barrier, as 
computed far from the origin where the barrier is sup- 
pressed, is — j7r/2. Therefore, we write the phase near 
the boundary r ~ i? as 



$(r) 
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dr' -j7:/2, 



(7) 



which is valid when the classical turning point <^ R. 
Furthermore, at the boundary we have $'(_R) = 



se 



eR/{R- 



The latter is computed by im- 



posing continuity of the metric at the boundary (where 
all but the energy term in Eqns. dU and ([5]) are sup- 
pressed). Defining 



A(r, 



(8) 



the logarithmic derivative of the interior wavefunction at 
r — R can be written as 



0'(r) 



eR 1 
cot [eA(r,) - j7r/2] 



R-r 



R 



(9) 



In the black-hole limit — 0, therefore from Eqn. ([8|) we 
see that as — >■ i?, A(rs) tends to infinity. In the quan- 
tum case this is because A(rs) is related (but not equal) 
to the total phase accumulated by the particle on the 
interior. This phase is large due to the wave function os- 
cillating many times on the interior as the particle moves 
rapidly in the strong field. However this integral also 
gives the classical time that a massless particle [ds^ = 0) 
spends on the interior, which goes to infinity in the black 
hole limit. 

It is worth noting that the phase in Eqn. ([7]) evalu- 
ated at the boundary does not correspond to the total 
phase accumulated over the interior as this requires ex- 
plicit knowledge of rt- 

However, the discrepancy with the total phase is only 
a constant offset in the limit — >■ 1, and (O gives rise to 
a good approximation to the wavefunction and its deriva- 
tive at the boundary. 



IV. GENERAL EXTERIOR SOLUTION 

The metric on the exterior of a static massive body is 
given by the standard Schwarzschild metric, which yields 
the following radial differential equation 



dr"^ 



1 



dr 



eV2 
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(r - rsf - ^s) - rs) 



(r) - 0. (10) 



A. Region I (r ~ rs) 

For e ^ (r—rg) (later verified in our region of interest), 
Eqn. pH)) becomes 



dr'^ 



1 

r j dr 

jU + 1) 



(r-Ts) r{r - Ts) 



0, 



(11) 
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which has the exact solution 



Hr)^a.[-) Pj-f^(l-^) (12) 
+ /3i (^)'^' 2F1 (j-s + l,j + ,s + l,2i + 2;^) , 



where P^'^\x) represent the Jacobi polynomials and 
2-^1(0, &, c; z) is the Gaussian hypergeometric function. 

B. Region II (r > r^) 

In this region we assume low energy in Eqn. (jlOp as 
before and make the substitutions 4'{r) — x(»')/r and 
r = p/e. This yields 



2er,, + 



P 



p- 



x = o. 



(13) 



The corresponding solution in terms of the regular and 
irregular Coulomb wave functions reads 

011 (rj = . (i4j 



By matching the coefficients of (|T2|) and (|T4|) in the limits 
r ^ Ts and er <^ 1 simultaneously, using the appropri- 
ate asymptotic forms of these solutions [13| , we arrive at 
the following relationship between the coefficients of the 
solutions in Regions I and II as 



a2 



P2 Pi {j - s)\{j + s)!C2(2j + l)(£r,)2j+i 
2je^'''''l'^\V{j + l-i£rs)\ 



(15) 



C 



(2j + 1)! 



S-MATRIX 



The solution to Eqn. dTOt at large distances can also 
be written in terms of outgoing and incoming waves as 



(16) 



where 



— er + er^ ln(2£r) + S^' — J7''/2 and 
= arg [r (j + 1 + ierg)] is the Coulomb phase shift. 
This allows us to write the scattering matrix as 



B, 



(17) 



Imposing the condition er ^ 1 in Eqn. p^ . we have 
the following asymptotic form of the wavefunction in re- 
gion II 



a2 sin{z) + P2 cos{z) 



(18) 



which, upon matching to (jl6p . gives the scattering matrix 
as 



S 



1 



ia2 

h_^2^5j 



1+ ^ 



(19) 



VI. MATCHING OF WAVEFUNCTIONS AT 
BOUNDARY AND RESONANCE ENERGIES 

Matching the logarithmic derivatives of the exte- 
rior (jl2p and interior Q wavefunctions at r = i? and 
taking the black hole limit r<; — > i? gives 



ai {~iy-^T{2j + 2) tan (eA(r,) - j7r/2) 



Pi 



sR r(j - s + i)r(j + S + 1) 



(20) 



Resonances occur at energies where the absorption cross- 
section is maximized, i.e. S* = — 1. This is achieved in 
Eqn. (HH) by setting a2l P2 = 0, which by Eqn. ([15]) is 
equivalent to ai/Pi — 0. Setting ai/Pi to zero in ((20|) 
results in the resonance condition for the energy 



"7^ + J2 

A(r.) ' 



(21) 



where n = 1,2,.... Note that these resonance energies 
strongly depend on the value of j, but differ marginally 
between s — Q and s = 1. 



VII. RESONANCE WIDTHS 

The full resonance condition corresponds to complex 
energies e = e„ — ir„/2 for which the scattering matrix 
([T9)) has a complex pole. This occurs when 



1 



ia2 



0. 



(22) 



On resonance ai/Pi — allowing us to expand this ex- 
pression near e = £„ as 



Pi 



is) 



ai 




d 




Pi 









(23) 



At the complex pole of the scattering matrix, 
£ — En — i^n/'^, this gives 



Pi 



d_ 



-iVr 



(24) 



where we have only substituted for energy where it ap- 
pears explicitly in Eqn. Using ([15]), ^ and 
(1^^ (with e e„, Ts i?), the resonance widths are 
given by 



r. 



(25) 



where we have assumed eh-{rg) 3> i7r/2, i.e. large n (see 
(1211) ) , which is already assumed when invoking the semi- 
classical approximation. 
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VIII. ABSORPTION CROSS-SECTION 

As discussed in Section Hill in the black hole limit we 
find that A{rs) tends to infinity. Thus by Eqns. (I?!]) and 
(PS)) both £„ and r„ tend to zero in the limit — i? 
for any fixed, finite values of n and j. However, the 
ratio r„/£) remains constant, where D — Sn+i — £n is 
the spacing between adjacent levels. This allows us to 
use the optical-model (energy-averaged absorption cross- 
section) This is obtained by averaging over a small 
energy interval containing many resonances and reads 



-opt 



E 

j-s=0 



D 



■(2j + l) 



(26) 



Substituting and (l25|) into ([Ml), gives 



E 



47r(2j + l)C^e^m^^+^{{j - sy.fjjj + 5)!)^ 



(27) 

which is independent of A(rs) and thus of the interior 
metric. 

By setting s = or s = 1 in Eqn. (j27p we recover 
the low-energy, total absorption cross-sections for mass- 
less scalar particles [1] or photons 9, 10], respectively, 
incident on a Schwarzschild black hole. 



IX. SPECIFIC INTERIOR METRICS 

In this section we present calculations involving two 
specific interior metrics that allow the — )■ i? limit 
to be taken: the Florides \XM ^^.^ Soffel [ll] metrics. 
Specifically, we verify our analytic solutions with nu- 
merically calculated resonance widths and energies via 
the short range phase shift 6(e). To calculate S{e) we 
solve the second-order differential equation numeri- 
cally, for given e'^ and e^, with the boundary condition 
(j){r — 0) ~ using Mathematica [l^. This solution 
provides a real boundary condition for the exterior wave 
function a.t r = R. (We set i? = 1 in the numerical cal- 
culations). Equation (|10p is then integrated outwards to 
large distances r ^ rg- In this region Eq. (jlOl) takes the 
form of a non-relativistic Schrodinger equation for a par- 
ticle with momentum e and unit mass in the Coulomb 
potential with charge Z = —rgS^. Hence, we match the 
solution with the asymptotic form [l| 

X{r) cx sin[er - (Z/e) In 2er + Sc + S - 311/2] (28) 

and determine the short-range (numeric) phase shift 5. 
This phase possesses steps of height tt at the resonance 
positions e„. We fit the step profile of an individual res- 
onance to the Breit-Wigner function 



arctan 



(29) 



r„/2 

where 5„ is a constant, from which we extract the numeric 
resonance widths and positions r„ and e„. 



0.9980 0.9984 0.9989 



0.9993 



0.9998 



FIG. 1. Energies of the n = 2 resonance in the Florides metric 
with s = J = 1. Closed circles indicate numeric data, the solid 
line indicates analytic e„ given by Eqn. (|21|) with h-ifa) given 
by Eqn. ([3lJ- 



A. Florides Interior 



The Florides metric is characterized by 



3/2 



This in turn leads to 



i?3 



(30) 



A(r. 



V2F(V4)F(5/4)(l-r3/i?)3/4 
1.198 (l-r,/i?)-3/^ 



(31) 



The resulting resonance energies and widths are com- 
pared with their numeric counterparts in Figures ([1} and 
([2]) respectively. 



B. SofFel Interior 

The Soffel metric is characterized by [lH| 

r,(l-rVi?2) 



exp 



2i?(l - Ts/R) 



(32) 



with e'*'^''^ equal to that of the Florides case. This in turn 
leads to 



Ts/R 



4(1 -r,/i?) 



(33) 



Analytic and numeric e„ and F„ for the Soffel metric are 
compared in Figures ^ and Q respectively. 



0.9980 0.9984 0.9989 0.9993 0.9998 

FIG. 2. Widths of the n = 2 resonance in the Florides metric 
with s = j = 1. Closed circles indicate numeric data, the 
solid line indicates analytic r„ given by Eqn. (I25|l with A(rs) 
given by Eqn. pT|) . 



"'0.950 0.955 0.960 0.965 0.970 

FIG. 3. Energies of the n = 4 resonance in the Soffel metric 
with s = j = 1. Closed circles indicate numeric data, the 
solid line indicates analytic e„ given by Eqn. (|2ip with A(rs) 
given by Eqn. ((33)) . 



X. CONCLUSIONS 

The problems of the scattering of low-energy, massless 
spin-0 and spin-1 particles from a massive, static, spher- 
ical body have been considered. We have shown that 
such scattering is characterized by a dense set of long 
lived resonances. Capture to these long-lived states gives 
rise to effective absorption in a purely potential scatter- 
ing setting. In the black hole limit the cross-section for 
absorption exactly equals the cross-section in the pure 
black hole case (for low energy) . Thus scattering of pho- 
tons (and massless scalars) by a near-black-hole object 
resembles black hole absorption. 

We thank G. F. Gribakin for useful discussions. This 
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FIG. 4. Widths of the n = 4 resonance in the Soffel metric 
with s = j = 1. Closed circles indicate numeric data, the 
solid line indicates analytic Fn given by Eqn. (|25|l with A{rs) 
given by Eqn. (f33)) . 
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